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ABSTRACT 


The Buffalo District, Corps of Engineers’ Lake Erie Wastewater Management 
Study and Heidelberg College's Water Quality Laboratory supported a tributary 
water quality monitoring program from 1974 to 1980 of the major United States 
tributaries to Lake Erie. For selected streams, sampling stations were 
established at the furthest downstream U. S. Geological Survey flow gaging 
station. This program was designed to monitor concentration changes 
occurring in association with increased streamflow with primary emphasis 
placed on sampling during high flow conditions (when the bulk of mass 
transport occurs). This report presents the authors' best estimate of the 
loadings for several nutrients to Lake Erie using collected data. A calcula- 
tional procedure employing stepwise concentration functions for fixed flow 
intervals was developed. Intervals which best represented the full range of 
streamflows were selected and used in conjunction with the daily flow 
records and the concentration functions to calculate monthly loadings. 
Loadings for unmeasured areas were extrapolated from the loadings for adja- 
cent areas. The loadings were then combined with municipal and industrial 
point source loading information and Lake Huron outflow data to produce a 
point, nonpoint, and total nutrient load summary for the Detroit~-St. Clair 
Basin and the three basins of Lake Erie on a monthly and yearly basis for the 
period 1970-1980. The nutrient loadings presented include total phosphorus, 
orthophosphate, nitrite-nitrate, ammonia, organic nitrogen, chloride, and 
silica. A three-basin, in-lake, total phosphorus model is presented and in- 
lake concentrations are modeled using the calculated 1970-1980 loadings. The 


results are compared with measured in~-lake values. 
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INTRODUCTION 
Public Law 92-500 authorized the U. S. Army Corps of Engineers to develop a 
wastewater management program for Lake Erie. The study identified the devel- 
opment of reliable nutrient loads to the lake as a major task. During the 
initial stage, a close examination of existing water quality records for the 
major Lake Erie tributaries was made. A program of bimonthly sampling had 
been carried out on the south shore tributaries from 1967 to 1969 by the 
Federal Water Pollution Control Agency (FWPCA 1969). State and previncial 
agencies carried out monthy sampling programs on some eributartes aiid the 
U. S. Geological Survey (USGS) also sampled some tributaries monthly. These 
programs with the exception of FWPCA's, were designed to monitor ambient 
water quality and not the mass transport of nutrients to the lake. The FWPCA 
report concluded, “Widely variable flow conditions were encountered at the 
time of sampling. Generally, concentrations did not show consistent changes 
with changes in flow; however, major loadings did occur during high 
discharges . The limited frequency of sample collection, particularly during 


periods of high flow, could bias the data.” 


Many agencies report total phosphorus loads to Lake Erie based upon the 
monthly sampling. Rarely are the loadings for other parameters reported. 
Total phosphorus loads to Lake Erie are reported annually by the Water 
Quality Board of the International Joint Commission. Each political 
jurisdiction reports tributary, municipal, and industrial loads. However, 
the sampling programs for the tributary loading estimates are often inade- 
quate. Usually 12 samples per year are collected, sometimes more, but 
usually less. Sampling programs at these frequencies are likely to 


underestimate the annual loads because most of the phosphorus transport 


occurs within a few days of high flow and these events are not usually 
sampled (Yaksich and Verhoff 1983). Since 1980, the Great Lakes National 
Program Office has supported event sampling programs which supplement the 
local jurisdictional data. Sampling frequencies have recently increased for 


some Canadian tributaries. 


An event sampling program was established by the Lake Erie Wastewater 
Management Study in the fall of 1974. This sampling program documented 
established tributary loading to the lake and the data were used to define 
nutrient loads to Lake Erie as presented within this report. These loadings 
are the authors’ best estimate for several nutrient parameters. These para- 
meters were chosen primarily because they are important to the lake's welfare 
and also because they were needed to interface with lake modeling efforts 
carried out at Manhattan College under the sponsorship of the USEPA's Large 


Lakes Laboratory at Grosse Isle, Michigan. 
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SAMPLING STRATEGIES AND ASSOCIATED ERRORS FOR LOADING ESTIMATES 


For event response rivers, such as those draining into Lake Erie, most of the 
nutrient and material transport occurs during storm events and the spring 
runoff period. Some materials such as total phosphorus and suspended sedi- 
ment increase in concentration with increasing flow, while others such as 
chloride tend to decrease with increasing flow. Other nutrients such as the 
nitrogen species show no definite relationship to changes in flow, although 
they may show some increase for increasing flows. Even where nutrient con- 
centrations decrease with increasing flow, the nutrient loading increases 
since the relative increase in flow is much greater than the relative 
decrease in concentration. Figure |] shows the resultant total phosphorus 
transport for a storm event. It is during the high flow events that the bulk 
of the material transport to Lake Erie takes place. Hence, it is imperative 
to have good sampling during this period in order to accurately define the 
nutrient transport. 


Fixed Interval Sampling 


The monthly sampling programs carried out by most agencies in the U. S. are 
examples of fixed interval sampling. In these programs, water quality 
samples are taken at predetermined times, at intervals such as once or twice 
a month. Treunert et al. (1974) examined the effect of sampling frequencies 
on the determination of the annual phosphorus load of small to moderate 
streams (Wahnbach River in Germany) with low flows and no relationship 
between total phosphorus concentration and flow (stable response river). 
They examined sampling frequencies of 1 to 29-day intervals between measure- 


ments. They also examined four different averaging methods to determine 
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Figure 1 - Daily Total Phosphorus Loads in the Maumee River at Waterville, 
1975. 
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annual loads. They concluded that use of continuous flow measurement by 
means of a level recorder leads to annual total phosphorus loads two to four 
times more accurate than discontinuous flow measurements with a simple 
measuring weir. In addition, they concluded that sample intervals of 14 to 
21 days (17 to 26 samples/year) and continuous flow measurements with a level 
recorder are required to keep deviations within 20 percent of the reference 


value. 


In contrast to the above results, Hetling, et al. (1976) concluded that a 
fixed interval sampling is of little value to obtain annual loads unless 
daily or every other day samples are taken. With sample intervals of 14 to 
21 days (17 to 26 samples/year) deviations of 700 percent were obtained for 
particulate phosphorus and 1,700 percent for suspended sediment for an event 
response river. They further concluded that a sampling strategy based on 
sampling high flow events was necessary to obtain reliable annual loads for 


an event response river. 
Event Sampling 


Yaksich and Verhoff (1983) examined high flow sampling strategies for event 
response rivers in the Lake Erie Basin (Maumee, Portage, Sandusky, and Huron 
Rivers). They examined total phosphorus, soluble reactive phosphorus, and 
sediment transport. They recommended a sampling program that included two to 
three of the largest events with 15 to 20 grab samples collected over each 
hydrograph, and an additional 5 to 10 samples collected during steady flow. 
They concluded that this program would yield a load estimate with an error of 
10 to 20 percent. They also determined that if high flows were not sampled 


loads would be underestimated by 15 to 30 percent. 


LAKE ERIE WASTEWATER MANAGEMENT STUDY WATER QUALITY SAMPLING PROGRAM 


A review of current research into river system chemical dynamics indicated 
that a significant new understanding of mass transport of pollutants had been 
developed. Two major studies, one in the Sandusky River Basin of Ohio (Baker 
and Kramer, 1973) and the other in the Brandywine River Basin of Pennsylvania 
(Cahill, et al. 1974) had utilized a sampling program of synchronized flow 
and chemical measurement at selected river sampling stations in an effort to 
develop accurate mass transport over a period of a year. By concentrating on 
the measurement of chemistry during storm runoff, these studies indicated 
that a significant portion of the total chemical pollutant transport out of a 
river system occurred during periods of high flow. More importantly, they 
showed very distinct increases in concentrations of certain water quality 
parameters during high flow, contrary to the dilutional effect that had been 


observed over the years with the standard sampling programs. 


Based on this information and because of the inadequacy of the available che- 
mical record for Lake Erie tributaries, a sampling program was designed by 
the Lake Erie Study staff during the autumn of 1974 to measure flow and che- 
mistry on selected tributaries. The initial stream selection was based pri- 
marily on the quality of flow records available, since the availability of a 
continuous flow recorder station at the sampling site would facilitate the 
measurement of stream flow changes with corresponding changes in water che- 
mistry. The longer the flow record, the more accurate the flow “rating 


curve.” 


It should be noted that this sampling program was designated on the basis of 


anticipating increased concentrations of certain pollutants, such as total 
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phosphorus, occurring in association with increased stream flows. For this 
reason, primary emphasis was placed on sampling high flow conditions, since 
the occurrence of increased concentrations during high flows would substan- 
tially increase mass transport. Historical chemical records, when 
available, were confined almost entirely to middle and lower flow regimes 


with virtually no data for the high flow range. 


Ten sampling stations were selected for monitoring during 1975. They 
measure runoff from 10,423 square miles or 51 percent of the total Lake Erie 
drainage area in the United States. Heidelberg College's Water Quality 
Laboratory, the city of Cleveland's Water Quality Program, and the State 
University College of New York at Fredonia were contracted for sample collec- 
tion and analysis during 1975. The Toledo Metropolitan Area Council of 
Governments contracted with Heidelberg College to continue stations on the 


Maumee, Portage, Sandusky, and Huron Rivers from June 1975 to June 1976. 


The Michigan Department of Natural Resources, Heidelberg College, and the 
Great Lakes Laboratory were contracted by the Corps of Engineers for 1977. 
The U. S. Geological Survey (USGS) also collected some samples for analysis 
by Heidelberg College. EPA's Athens Laboratory and the Rockefeller 
Foundation provided support on the Maumee, Portage, Sandusky, and Huron 
Rivers for the second half of 1977 and 1978. The Great Lakes National 
Program Office contracted Heidelberg College and the USGS for sample analysis 
and collection in 1978 and 1979. Table 1 shows the tributaries that were 
sampled and used for the loading calculations. Many more sampling stations 
were established than those used here, mostly upstream of these stations 


(Corps of Engineers, 1982). For the Canadian side, water quality data were 
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obtained through the Ontario Ministry of the Environment and flow data 


through Canada Water Survey. 


Samples collected were analyzed for the following parameters: total and 
orthophosphorus, nitrite~nitrate nitrogen, ammonia nitrogen, chlorides, 
dissolved silica, suspended solids, and conductivity. Twenty percent of the 
samples were analyzed for total kjeldahl nitrogen, inorganic nitrogen, and 
iron. Analytical methods are summarized in “Lake Erie Basin Tributary Mouth, 
Water Quality Data" (Corps of Engineers, 1977). The sampling procedures and 
results for the northwestern Ohio River have been described in detail by 
Baker (1982). Table 1 shows the sampling dates and number of samples taken 
on the U.S. streams under direct support from the Corps of Engineers. The 
map codes correspond to locations in Figure 2. Data for these stations were 
augmented by data collected and analyzed through USEPA support and data from 
the USGS. The total number of samples used in the loading calculations is 
also shown in Table 1. Stream flow data were provided by the Michigan, Ohio, 


and New York Divisions of the U. S. Geological Survey. 


LOADING CALCULATIONS 


Tributary Loads 


For nutrient loading calculation purposes, the Lake Erie drainage basin was 
broken into 80 subbasings. These subbasins correspond to the major tribu- 
tarles flowing into Lake Erie, the Detroit River, Lake St. Clair and the St. 
Clair River, and the adjacent areas between these tributaries with direct 


drainage to the lake. These 80 subbasins are shown in Figure 2. 
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For the 19 sampled main stem streams on the U. S. side of Lake Erie and the 


12 Canadian streams, nutrient transport loadings were calculated on a monthly 
basis for the period 1970-1980. The method used is called the flow interval 
method and has been described by Verhoff, et al. (1980). The flow interval 
method divides the expected discharge range into groups. Within each group, 
a flow weighted mean concentration is calculated. Daily flows are used with 
the concentration which was calculated for the interval in which the flow 
belongs to calculate the load for the day. Daily loadings were then sum- 


marized for monthly loads. 


Figure 3 shows a plot of concentration versus flow. The concentration may 
increase, decrease, or not be affected by increasing flow. One method of 
deriving a concentration-flow relationship is to approximate this relation by 
a least squares curve as shown in Figure 3. As can be seen, this does not 
necessarily define the relationship very well, in some cases. However, a 
much better approximation can be achieved when the flow is broken up into 
convenient groups and flow weighted mean concentrations are calculated. The 
flow groups may be picked by any convenient method such as taking the maximum 
flow and dividing it into a certain number of intervals. For the calcula~ 
tions and results presented in this paper, a unit area flow was chosen. This 
was done for two reasons. First of all, the intervals would be uniform for 
all years and would not be related to any one year's high flow, and there- 
fore, results for different years could be compared. This comparison allows 
known data to be substituted in years which had missing data. Secondly, the 
flow weighted mean concentrations can be compared for different streams. The 


unit area flow intervals chosen were not uniform in size. They were small on 
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the low end and large on the high side, and as many as 30 intervals were 


used. 


If enough data points existed for a year and they were high flow oriented, 
the flow weighted mean concentrations were calculated for each flow interval. 
If a given year lacked enough data to be used in the calculation, these data 
were combined with data from an adjacent year and used in the calculation. 

In the case in which data were totally lacking for a year, flow weighted mean 
concentrations averaged over the years for which there were data were used. 
This was possible because the unit area flow groups were the same size for 
every year. If only one or no points existed for a given interval, the flow 
weighted mean concentration for that interval was approximated from two adja- 
cent intervals in that year. At least two data points per interval were 
required. In the situations where high flow data were missing the con- 
centration from an adjacent year for the appropriate unit area flow interval 
were used. In some cases, not enough samples were taken or no samples were 
taken during high flow for the entire period of record. In these cases a 
least squares regression of total phosphorus with discharge that was genera- 
lized for the Lake Erie watershed was used to define the daily phosphorus 


loadings (Verhoff, et al., 1980) 


It should be noted that large seasonal differences in concentration versus 
flow relationships are present and well documented for the intensively 
Studied basins. Since correspondingly detailed data were not available for 
most monitored streams, annual rather than seasonal concentrations versus 
flow relationships were used. This would tend to overestimate winter and 


early spring loading and underestimate late spring and summer loading, but 
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these should average out for the year. Using this method monthly loads were 
calculated for total phosphorus, soluble reactive phosphorus, nitrate-nitrate 
nitrogen, ammonia nitrogen, organic nitrogen, suspended solids, chlorides, 


and soluble silica. 


Point Source Loading Calculations 


Point source nutrient loadings for the 80 defined basins were calculated 
from the computer summary of Volume II, Appendix A of the Lake Erie Wastewater 
Management Study Preliminary Feasibility Report (1975) updated thru 1980. 
These data were originally compiled from local, regional, and Federal data 
bases for both the U. S. and Canadian sides of the lake. Most of the infor- 
mation was obtained from EPA, State, and local permit records. It was con- 
tinuously updated through yearly reports, especially those from the Water 
Quality Board of IJC, when these became available. The information from each 
facility in the 80 identified watersheds was compiled with average yearly 


flows and accompanying concentration data. 


Where past data were not known or unavailable, the most recent data were 
assumed to have remained constant over the preceeding years. An exception to 
this assumption is treatment plants known to have implemented phosphorus 
removal. Total phosphorus concentrations were assumed to be 6mg/l before 
phosphorus removal. When no chemistry data were available, as in the case of 
small package plants, a concentration was assumed for that plant over all 
years. Table 2 lists the seven parameters of interest with default concen- 
trations which were assumed when no data were available. In most cases, 


wastewater treatment plants with flows greater than one million gallons per 
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day reported values for at least total phosphorus. The default con- 
centrations were based on the discharge from other municipal plants in the 
Lake Erie Basin. For many of the very small facilities, there was a lack of 
both flow and chemistry data. Because of the insignificant contributions 
these small facilities had on Lake Erie as a whole, they were eliminated from 
the point source file. Nutrient loadings for the seven parameters are calcu- 
lated for each facility and summed for each basin above and below the flow 
gaging station. Calculations for the point source loads were made using the 
average concentration for a year times the average flow for that year. The 
monthly point source loadings were obtained by proportioning the yearly loads 
by the number of days in the month. 


Table 2 - Default Concentration Values for Effluent Concentrations 
from Wastewater Treatment Plants 


Total Phosphorus 6-0 mg/l as P 


Orthophosphate 3.0 mg/l as P 
Nitrate-Nitrite Nitrogen 6.6 mg/l as N 


Ammonia Nitrogen mg/l as N 


Organic Nitrogen 2.33 mg/l as N 
Chloride 160 mg/l as Cl 


Silica 5 mg/l as Si02 
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Lake Huron Outflow Loadings 


Lake Huron input to the St. Clair River was handled as a separate entity 
because of the nature of the water quality data that were available. Data 


were obtained from the Ontario Ministry of the Environment in the form of 
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panel measurements across the river. For each panel, flow and concentrations 
were reported which represented average values for the given panel. From 
these, a total flow and a panel flow weighted mean for the parameters of 
interest were calculated for each date information existed. From the yearly 
records, flow weighted mean concentrations for the year were calculated. In 
some cases, the yearly record for certain parameters was spotty. If it was 
judged there were not enough values to have a meaningful flow weighted mean 
concentration, then data for that parameter were combined with all other 


year's data for that parameter and these were used for the caJculations. 


The average daily loading for the head of the St. Clair River was obtained by 
multiplying the flow weighted mean yearly concentration by the reported 
average daily flow per month. This daily loading was prorated by the number 


of days in the month to yield a monthly loading. 


Point/Nonpoint Nutrient Load Summary 


At this point, the St. Clair River loadings, the tributary loading cal-sula- 
tions, and the point source calculations must be brought together to produce 
a point/nonpoint load summary for each of the 80 designated basins on a 
monthly basis. Nonpoint or diffuse loadings must be determined for the 80 


regions both above and below the gage where they exist. 


For the 3] tributary streams, the monthly upstream point sources above the 
gage are subtracted from the tributary loads. This defines the above gage 
nonpoint loading. It should be noted, in the case of indirect point sources, 
during low flow months much of the phosphorus input would be tied up in the 


stream sediment during these months. In some cases the tributary flow is so 
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low that a negative nonpoint source yield is calculated. In this case the 
nonpoint yield is taken as zero and point source load is set equal to the 
tributary load. Dividing the nonpoint source load by the drainage area above 
the gage defines the above gage unit area nonpoint loading. The downstream 
nonpoint load below the gage is calculated by multiplying the unit area load 
by the area below the gage. Loadings for the ungaged areas between the gaged 
streams were calculated by taking the average unit nonpoint area loads from 
the two adjoining gaged streams, and mitiplying that average by the ungaged 
area. This same extrapolation procedure was used for streams that did not 
have a water quality sampling program. The total loading for each of the 80 
basins is obtained by summing point and nonpoint loadings above and below the 


gages. 


Figure 4 is used to illustrate the extrapolation of nonpoint loads for 
measured basins to those basins not measured. Through the water quality 
sampling program, the tributary loads for areas B and D are defined. 

Nonpoint loads for these areas can be calculated by subtracting off the point 
source loads. Nonpoint loads for the area below the sampling points, A and 
C, can be defined in terms of the unit area loads for areas B and D, respec- 
tively. For example, the nonpoint loads for A equals the unit area load for 
B multiplied by the area of A. The nonpoint load for E can be defined as the 
area of E multiplied by the average of the unit area loads for B and D. 


Similarly, loads for F and G can be defined. 


To produce monthly and yearly loading for the three basins of Lake Erie and 
the Detroit/St. Clair River system each of the 80 basin values are added 


together for the appropriate class. Basins 1-10 and 75-80 are in the 
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Table 3 - Yearly Summary for 1970 - Basin Loading in Metric Tons per Year 


vo se oe LEI Oe eo ee ee ee 
: St. Clair : Western : Central : Eastern : 
: Basin : Basin : Basin : Basin : Total 
Total Phosphorus : : : : : 
Lake Huron : 2,502 ; : : : 2,502 
U.S. Nonpoint : 307 : 2,389 :; 1,244 ; 296 : 4,236 
Canadian Nonpoint 966 27 : 210 : 608 1,811 
Total Nonpoint : 1,273 : 2,416 : 1,454 ; 904 : 6,047 
U.S. Point 6,622 1,511 3,035 : 538 11,706 
Canadian Point : 91 ; 0 ; 12 : 90 : 193 
Total Point : 6,713 : 1,511 : 3,047 628 ; 11,899 
U.S. Total ae 7 ee ea ee 
Canadian Total : 1,057 : 27 : 222 ; 698 ; 2,004 
Basin Total : 10,488 3,927 : 4,501 1,532 20,448 
Ortho Phosphorus : : : ; 
Lake Huron : 1,129 ; : : : 1,129 
U.S. Nonpoint : 86 : 309 : 123 : 1 519 
Canadian Nonpoint ; 244 38 108 : 116 : 506 
Total Nonpoint : 330 347 : 231 : 117 : 1,025 
U.S. Point : 3,319 : 665 : 1,341 261 : 5,586 
Canadian Point : 47 0 : 10 : 47 104 
Total Point : 3,366 : 665 : 1,351 : 308 : 5,690 
U.S. Total : 3,405 . 974 : 1,464 : 262 : 6,105 
Canadian Total : 291 38 : 118 163 : 610 
Basin Total : 4,825 : 1,012 : 1,582 : 425 : 7,844 
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Figure 5 Monthly Total Phosphorus Loads to Lake 
Erie, 1970-1980. 
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1980 


drainage area for the Detroit/St. Clair Basin. The total input from this 
basin to the Western Basin of Lake Erie is the sum of these plus the Lake 
Huron input. The total added tributary loading to the Western Basin, 
excluding the Detroit River input, is the sum of basins 11-20 and 72-74. The 
total loading to the Western Basin is the sum of these and the Detroit River 
input. The tributary loading to the Central Basin is the sum of basins 21-45 
and 62-71, while that for the Eastern Basin is the sum of 46-61. For any 
modeling purposes, it must be remembered that these are the external sources 
to the lake and the Central and Eastern Basins would also receive inputs from 
these upstream basins. For these calculations, atmospheric loadings and any 


phosphorus .added due to shore erosion were not added. 


Table 3 illustrates part of the yearly summary for 1970. A complete set of 
the yearly and monthly loading summary can be found in Appendix 1 and 2 of 


this report. 


RESULTS AND DISCUSSION 


Monthly loadings to Lake Erie exhibit considerable variation for the ll-year 
period, 1970 to 1980. An example of this variation is the total phosphorus 
loads shown in Figure 5. In general, large increases in loading are seen to 
occur during the spring runoff. Loads decrease during the summer and usually 
remain low during the fall. Despite the large variation in total phosphorus 
loads that occurred, a decrease in the total phosphorus for 1970 to 1980 can 
be seen from Figure 5. The decrease in the minimum monthly loads reflects 


the reductions in point source loading. 
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Figure 6 shows a comparison of the monthly loading between years 1977 and 


1978 for the parameters total phosphorus, soluble orthophosphate, nitrite- 
nitrate nitrogen, ammonia nitrogen, organic nitrogen, chloride, and silica. 
Notice the general spring runoff peaks during March and April with smaller 
peaks through the year representing storms during those months. 1977 was 
illustrated because of the large storm that occurred in December of that 
year. It should be noted that parameters whose concentrations decrease with 
increasing flow (soluble orthophosphate and chloride) still experience the 
maximum loading during the runoff period. However, they show about a twofold 
increase while total phosphorus shows a fourfold increase during high flow 


periods. 


Figure 7 shows the yearly loading trends for all the parameters for the 
years 1970-1980. The curve labled Detroit River is the loading at the mouth 
of the Detroit River. it represents the sources added to the Detroit/St. 
Clair Basin along the way from Lake Huron to the mouth of the Detroit River, 
plus the Lake Huron load. The Western Basin load is the sum of all sources 
to the Western Basin excluding the Detroit River. The Central Basin ioad 
does not include the input from the Western Basin and the Eastern Basin load 


does not include the input from the Central Basin. 


The total phosphorus loading to the lake as a whole starts to decline sharply 
in from 1973 to 1980. This is due mainly to the sharp decrease in the 
sources to the Detroit Basin as can be seen in Figure 7a. The same general 
trend is exhibited in Figure 7b for orthophosphorus. The nitrogen species 
(Figures 7c, d, and e), show more variability with somewhat erratic behavior. 


This is especially true with the nitrogen loading components into the Western 
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Figure 8 Total Phosphorus Loads to Lake Erie, 1970-1980; 
Point and Non-Point (Does not include shoreline 
erosion or atmospheric input.) 


Basin of the lake including Lake Huron. It is interesting to note, however, 
that the values for nitrite-nitrate, ammonia, and organic nitrogen loadings 
remains relatively constant for the Eastern Basin of Lake Erie. Overall, the 
nitrogen loading into the lake has increased over the period. The chloride 
load to the lake has decreased from 1970 to 1980 due to the drop in the sour- 
ces added to the Detroit River and a general decrease after 1976 for the 
Central Basin sources. Silica is sporadic like the nitrogen species, but 
does show a decrease in loading from 1975 on. This is basically due to large 
drop in the output from Lake Huron. It should be noted that the data for the 
St. Clair River is sparse of many years. Silica does show a seasonal 
variation and changes in the Lake Huron loading could be due to inadequate 


sampling. 


Figure 8 shows the changes in the point and nonpoint total phosphorus 
loadings to the lake as a whole. Lake Huron output has remained relatively 
constant. Both point and nonpoint show a somewhat decreasing trend, point 
more so than nonpoint. Nonpoint is sensitive to rainfall and hydrology and 
can be expected to be more variable from year to year. In the last half of 
the 1970's there has been a general trend in reducing sewage treatment plant 


total phosphorus output. 
THREE-BASIN LAKE ERIE TOTAL PHOSPHORUS MODEL 


A three~basin total phosphorus model for Lake Erie, (Corps of Engineers 1975 
and Rumer 1977) was used in conjunction with the total phosphorus loadings 
calculated from 1970 to 1980 to predict water quality changes over that 
period. The model is basically a six-compartment model with a water and 


sediment compartment for each basin. Inflow into each basin is equal to the 
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basin outflow. The water component is assumed to be continuously stirred and 
has a uniform total phosphorus concentration in each compartment. Figure 9 
is a schematic showing the three basins and the processes taking place. 

M, Tepresents the tributary and Detroit River loading to the Western Basin, 


while My and M3 represent tributary loading to the other basins. C and C, 


represent the total phosphorus concentrations in the water and sediment. 


The material balance equation for phosphorus in the lake water can be written 


as: dc M K(C3-C)A K,CA CQ 


—=s-4 


-—_—-— (1) 

dt Vv Vv V Vv 
where V is the lake basin volume assumed constant in cubic meters and A is 
the sediment area in square meters. The quantity dC/dt represents the annual 
change in lake phosphorus concentration required to balance the phosphorus 
budget. If the basin is at equilibrium and the phosphorus loading is main- 
tained at a constant, the lake phosphorus concentration remains unchanged and 


dC/dt equals zero. 


The material balance equation for phosphorus in the sediment interstitial 
water is given by: 


de Wg ge ge 


8 
Where V, is the interstitial water volume. Again, for equilibrium con- 


ditions, the quantity dC,/q+ equals zero. 


The above set of coupled equations can be solved for C and Cs as a function 


of time measured in years. The accuracy of the solution in quantitative 


—_— —_ er 2 a 


terms depends on the adequacy of the model formulation, correctness of the 
values chosen for the three rate constants (K,, Kj, and K3), adequate infor- 
mation on the phosphorus loading rate (M), and data on the lake character- 


istics (V, V,, A, and Q). Initial values of C and C, must also be known. 


There is an equation for each concentration of the water and sediment com- 
partments in each of the three basins. Five processes take place for each 
basin; a material transport to the basin, a material transport from the 
basin, a transport of phosphorus to the sediment, a transport from the sedi- 
ment in the water and a permanent loss from the sediment. The last three 
processes are controlled by the reaction parameters K), Kj, and K3, respec~ 


tively, which must be determined from existing data. 


If equilibrium conditions between the total phosphorus concentration and the 
loading are assumed (dC/dt = 0), the following two equations are obtained 
which relate the three rate constants. 


M 
K)K3 2 as @ 
AC,. (3) 


Kj (1-K3)C 
Ka = coke (4) 
2 Cg, - 
In this case, C, and Cy. are the steady-state or equilibrium values for lake 
and interstitial water concentrations, respectively. With knowledge of any 
one of the rate constants, the other two are specified by the above (3) and 


(4). 
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Figure 10 Lake Erie Total Phosphorus Concentrations 


A - Model Prediction for 1970 Loadings and 
Concentrations 

B = Model Prediction for Average 1970-1975 
Loadings and Concentrations 

Data Points and Error Bars Showing + 1 

Standard Deviation from Actual Mean 

Concentrations Are Presented. 
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For illustrative purposes and to check calculated basin concentrations with 
field data, two sets of reaction parameters were calculaicd and used in this 
model. In both cases, the values for Ky (transfer of phosphorus from 
sediment) were the same and these values were approximated from Project Hydro 
(Burns and Ross 1972) data as 1.0, 1.4, 1.0 m/yr for each of the basins. 
Using 1970 loadings and concentrations as equilibrium conditions, Case A 
parameters K, and Kj Were calculated as 59.1, 44.66, 30.63 m/yr and 0.872, 
0.594, 0.620, respectively. Similarly for Case B where 1970-1975 average 
loadings and concentrations were used as equilibrium conditions, the parame- 
ters are 79.87, 45.23, 13.17 m/yr and 0.905, 0.565, 0.385. These two cases 
were then used with the 1970-1980 loadings to predict in-lake total 


phosphorus concentrations for each of the basins. 


Figure 10 shows the model results plotted against measured in-lake total 
phosphorus concentrations. The plotted values with the bars are the ranges 
of the measured in-lake concentration which were reported by the Water 
Quality Board to IJC (1981) as measured by various agencies. As shown, the 
two cases predicted match the measured values fairly well. Only for the 
eastern basin after 1976 for Case B are the predicted concentrations outside 


the range for the measured data. 


The lake loadings presented in the report (total phosphorus loads and the 
others) were also used to calibrate, verify, and predict lake concentrations 
for all the parameters using a much more sophisticated model. This model was 
developed for the USEPA at Manhattan College. Calibration and verification 


of the model was very good and the predicted concentrations agreed reasonably 


well with measured data (DiTore and Connolly). 


CONCLUSIONS 
Total phosphorus loads to Lake Erie have decreased during the period 
1970-1980 as a result of phosphorus removal at wastewater treatment plants. 
The effect of these reductions can be seen in the lake concentrations and 
were predicted by a three-basin phosphorus budget model developed in the 
early 1970's. Soluble orthophosphate loads, chloride loads and silica loads 


have decreased over the period. Nitrogen species were highly variable and 


increased over the period. 
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Appendix 1 
Yearly Nutrient Loads to Lake Erie 
1970-1980 
‘ 
37 
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This appendix presents a yearly summary for calculated 
nutrient loads to Lake Erie for the years 1970-1980. 
The presented loads are tributary loads only and do 

Not account for shoreline erosion, atmospheric inputs, 
and sediment inputs. The loads are for total phosphorus 
as P, nitrite-nitrate nitrogen as N, ammonia nitrogen 
as N, organic nitrogen as N, chloride as Cl, and silica 
as SiO and are broken down into components for the 
Detroit/St. Clair River Basin, the Western Basin, the 
Central Basin, and the Eastern Basin of Lake Erie. 
These are further divided into U.S. and Canadian point 
and non point, total point and non point, and total 
U.S. and Canadian loads. 
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an ace reeset ae eee ey ES 


This appendix presents a yearly summary for calculated 
nutrient loads to Lake Erie for the years 1970-1980. 
The presented loads are tributary loads only and do 

hot account for shoreline erosion, atmospheric inputs, 
and sediment inputs. The loads are for total phosphorus 
as P, nitrite-nitrate nitrogen as N, ammonia nitrogen 
as N, organic nitrogen as N, chloride as Cl, and silica 
as SiO and are broken down into components for the 
Detroit/St. Clair River Basin, the Western Basin, the 
Central Basin, and the Eastern Basin of Lake Erie. 
These are divided into point and non point and total 
point and non point loads. 
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